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ABSTRACT 

We discuss the inclusive cross section to produce two minijets 
with a large separation in rapidity in high energy hadronic col- 
lisions. The contribution to the inclusive cross section from 
the exchange of a BFKL Pomeron is compared with the con- 
tribution from the exchange of two BFKL Pomerons, which is 
induced by the unitarization of the semi-hard interaction. The 
effect of the multiple exchange is studied both as a function of 
the azimuthal correlation and as a function of the transverse 
momentum of the observed minijets. 







1. Introduction 

One of the main topics in perturbative QCD is presently represented by semi- 
hard hadronic interactions, namely by hadron interactions with momentum trans- 
fer constant with energy but large enough to apply perturbation theory. One 
of the characteristic features of this kinematical regime is the large size of the 
corresponding cross sections, which, although in the perturbative domain, rise 
rapidly with energy. In fact, already at the energies of present hadron colliders, 
one may easily obtain semi-hard cross sections whose size is comparable to the 
total hadronic cross section[l,2]. At the partonic level, in a typical interaction 
configuration, one of the two interacting partons has a finite fraction of the par- 
ent's hadron momentum while the other one has a momentum fraction close to 
zero. The separation in rapidity of the two partons is therefore increasingly large 
with energy and, in the parton-parton cm. system, the transverse momentum 
exchange is small with respect to the longitudinal momenta. The Regge limit is 
then approached in semi-hard interactions not only in the whole hadron-hadron 
process but also in the underlying parton-parton interactions. 

When considering the large pt regime the momentum exchange is of the order 
of the incoming partons momenta. At the parton level such a large scale factor 
can be transferred only in a few interaction vertices and, as a result, the elastic 
two body parton collision is a good first order approximation to the elementary 
partonic interaction. In the semi-hard regime, since the semi-hard scale is small 
with respect to the total energy available, there are several parton vertices with 
momentum exchange of the order of the semi-hard scale. A consequence is that all 
semi-hard radiated gluons are to be taken explicitly into account for a proper fac- 
torization of the semi-hard component of the interaction. When the 2 — > n, rather 
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than the 2 — > 2, is the parton subprocess relevant to the semi- hard component 
of the hadronic interaction, a difficulty arises in constructing an inclusive cross 
section, where only few of the radiated partons are actually detected as mini-jets 
in the final state. In fact one is not allowed any more to use the lowest order 
tree diagram to represent the parton amplitude, since the tree level amplitude is 
singular in the soft and collinear limit. To avoid the infrared problem one faces 
when evaluating an inclusive cross section, one needs to keep virtual corrections 
explicitly into account and, as a consequence, the elementary subprocess acquires 
a non trivial structure. The problem has been addressed already several years ago 
in a series of papers by Lipatov and collaborators [3] . Lipatov's solution is the 
BFKL Pomeron: the partonic interaction is described by the exchange of a gluon 
ladder structure with vacuum quantum numbers in the t-channel. The s-channel 
discontinuity of the BFKL Pomeron represents the production of the semi-hard 
gluons. In the limit in which the transverse momenta are always negligible with 
respect to the longitudinal ones, the steps of the ladder are ordered in rapidity 
and dynamics is greatly simplified. Indeed the simplified kinematics lets one to 
isolate the two basic elements which build up the ladder: 
a- The gauge independent non-local vertices, which keep into account the domi- 
nant term, in the t/s — > limit, of the diagrams with gluon emission from all 
near-by lines, and 

b- the Reggeization of the t-channel gluons, which is the virtual correction that 
allows a solution to the infrared problem. 
The ladder structure can be iterated in the t-channel, which may be expressed as 
an integral equation, the Lipatov's equation. Lipatov's equation allows an analytic 
solution free from infrared (and ultraviolet) singularities. One obtains in this way 
an explicit expression for the cross section where two gluons interact producing 

2 



many gluons and two of them, the ones nearby in rapidity to the interacting 
partons, are observed. If y is the separation in rapidity of the interacting gluons 
and k a , kb are the transverse momenta of the observed ones, the inclusive cross 
section can be expressed as 



da 



d 2 k a d 2 kb 



f(k a ,k b ,y) 



k 2 



(1) 



where Ca = N c is the number of colors, a s is the strong coupling constant and 
f(k a , kb, y) is the inverse Laplace transform of the solution to Lipatov's equation. 
Actually: 



f{k a ,k h ,y) 
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where 4> is the azimuthal angle between the observed gluons, 



(2) 



LJ(v,n) = —2 M 



7T 



n I +1 



+ iv 



and 



dlnT(z) 

dz 



(3) 



(4) 



is the Digamma function. The inclusive cross section for production of two mini- 
jets, as a result of a BFKL Pomeron exchange, is obtained by folding Eq.(l) with 
the structure functions of the interacting hadrons A and B: 



da 



dx Adx Bd 2 k a d 2 k b 
where / e / / is the effective structure function 
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feff(xA,kl)f e ff(x B ,kl) 



da 



d 2 k a d 2 kb 



(5) 



f eff (x) = G(x) + [Qf(x) + Qf(x)] (6) 

/ 

namely the gluon structure function plus 4/9 of the quark and anti-quark structure 
functions with flavour /. The expression for the cross section in Eq.(l) correlates 
the azimuthal angle cj) with the distance in rapidity of the observed partons. The 
differential cross section in Eq.(l) may be easily integrated, at fixed, on k a and 
kb down to the lower cut off k m , which represents the threshold in transverse 
momentum that allows a parton to be observed as a minijet in the final state. 
This yields: 



rlfr (Cadi ) 2 1 ■ , r + °° p^O' 11 )?' 



which is a suitable expression to study the azimuthal correlation of the observed 
partons as a function of the rapidity difference y. Simple expressions may also be 
obtained for the cross section where the momentum of one of the two observed 
gluons has been integrated down to the lower limit k m : 



da _{C A a s f 1 f + °° ^ e^y / k a \ 2iv (g) 



d?k a 2n k^km J—oo 2iv + 1 \k m / 

and for the cross section where both observed gluons momenta have been inte- 
grated down to k m : 

The high energy behaviour of the integrated cross section is estimated by evalu- 
ating the asymptotic limit of the integral in Eq.(9) for large y[4\: 



7i(C A a s ) 2 exp[41n2N c a a y/ic\ 

& — ► — tzt^ T7T7 10 

2k l [7C(3)N c a s y/2} 1/2 
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where ( is the Riemann zeta function. Eq.(lO) shows a rapid growth, roughly linear 
with the parton-parton cm. energy, of the partonic cross section corresponding 
to the exchange of a BFKL Pomeron and gives a justification to the large size of 
the semi-hard cross section. 

The possibility to describe the elementary parton process by means of Lipa- 
tov's dynamics has been considered recently in a series of papers [5]. One of the 
main points of interest is the search for clear signatures of the underlying par- 
ton dynamics in the final state of high energy hadronic collisions. Correlations 
in transverse momentum and azimuthal angle, as a function of the distance in 
rapidity y of final state minijets, have been therefore estimated according to the 
expectations of the Lipatov's picture of the interaction as expressed by Eq.(7) and 
(8) [6]. On the other hand, to approach Lipatov's limit, one needs to keep the lower 
threshold of the transverse momenta of the observed minijets as small as possible, 
compatibly with the requirement of being still in the perturbative regime. As 
shown by Eq.(lO) smaller values of k m correspond to larger values of a. As a 
consequence of the larger probability of the elementary partonic intercourse one is 
therefore forced to take into account the possibility of having several elementary 
partonic collisions in each inelastic hadronic event, in order to implement unitarity. 

In the present paper, by assuming the validity of the AGK cutting rules in 
semi-hard interactions, we unitarize the semi-hard cross section and we derive the 
most general correction term to the inclusive cross section in Eq.(5). The correla- 
tions among the minijets observed in the final state are then estimated, considering 
the simplest possibility of multiple parton interaction, and are compared with the 
expectation from the single BFKL Pomeron exchange. The paper is organized as 
follows: in the next section the unitarity correction to the single scattering term is 
derived. In the following paragraph a few numerical estimates are presented, with 

5 



the purpose of giving an indication on the kinematical region where corrections 
can be expected to become sizeable. In the last section the general features of 
the unitarization of the single scattering term are summarized and a few general 
conclusions are drawn. 

2. General framework for muliparton interactions 

In order to approach the problem of multiparton interactions, with the pur- 
pose of obtaining for the inclusive cross section an expression which is more general 
with respect to Eq.(5), we find it appropriate to introduce a functional formalism, 
which keeps to a minimum level the occurrence of cumbersome expressions. As a 
preliminary step, we show how to derive all the relevant inclusive cross sections in 
the simple case of a single parton-parton collision with fixed fractional longitudinal 
momenta (x,x'). Let us introduce the functional 



where z is the argument of the functional and da n is the differential cross section 
to produce n partons with momenta (ki, . . . , k n ). Obviously the value of the 
functional for z = 1 is the semi-hard parton cross section a(x,x'). Actually 



All the inclusive cross sections are generated by taking an appropriate number of 
functional derivatives of the generating functional with respect to z[7]: 




(11) 



Q[x, x'; 1] = a(x, x') 



(12) 



da(x, x') %n 
dki . . . dk 



8z{ki) . ..Sz(k n ) 



5<d[x, x'; z] 



(13) 



6 



To obtain the inclusive cross section in the case of the actual hadronic collision 
a more elaborate analysis is needed. In the case of soft interactions multi-Reggeon 
exchanges are conveniently taken into account by making use of the AGK cutting 
rules [8]. Although no general proof of their validity is available in the case of 
semi-hard interactions, it has nevertheless been possible to show that the cutting 
rules hold for one of the components of the interaction which is leading in the 
large- s fixed-t limit [9]. If one assumes the validity of the cutting rules for semi- 
hard interactions, one is allowed to represent the semi-hard cross section oh as 
a probabilistic distribution of multiple semi-hard parton collisions[10]. The most 
general expression for ou requires however the introduction of the whole infinite 
set of multiparton distributions [11], which keep into account hadron fluctuations 
in the parton number: 

oh = J d 2 (3a H ((3) 

^)=/EE h w A } («!••• «») - Vir } («; - p . . . u' m - p) 

J ft. lib* 

n m 

n m (14) 

x{l- JJfJ[l «;■)]} JJdudu' 

i=l j=l 

Here the W^ k \ui . . .Uk) are the exclusive /c-body parton distribution, namely the 
probabilities to find a hadron in a fluctuation with k partons with coordinates 
u\ . . .Uk, ui = (bi,Xi) standing for the transverse partonic coordinate (pi) and lon- 
gitudinal fractional momentum (xi). (3 is the impact parameter between the two 
interacting hadrons and a(ui,u'j), represents the probability for the parton % of 
the A-hadron to have an hard interaction with the parton j of the S-hadron. The 
semi-hard cross section is constructed by summing over all possible partonic con- 
figurations of the two interacting hadrons (the sums over n and m) and, for each 
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configuration with n A-partons and m £?-partons, summing over all possible multi- 
ple partonic interactions. This last sum is constructed by asking for the probability 
of no interaction between the two configurations (actually Y\a=i Eljlitl ~~ )• 
The difference from one of the probability of no interaction gives the sum over all 
semi- hard interactions. cth(P) is then the probability to have at least one semihard 
parton interaction when the impact parameter in the hadronic collision is equal to 
(3. The semi- hard cross section is obtained by integrating the probability <jh((3) 
on the impact parameter. Analogously, the elementary semi-hard cross section 
a(x, x') is obtained by integrating the elementary interaction probability a(u, u') 
on the relative transverse coordinate b — b'. The expansion of on (P) as a sum on 
multiple interactions reads: 



^=/EE ^\ W A ] («!•■■ ««) ^ {<-(l...U> m - (3) 

n m 

xS Yl ( ja 1 ...a N (l-a N+1 )...(l-a Q ) 

IV— 1 V / 



S is a symmetrizing operator, which one may conveniently introduce taking ad- 
vantage of the symmetry of for permutations of the arguments[12], and the 
index iV counts the interactions which, for a given configuration with n A-partons 
and m B-partons, range in number from 1 to Q = nm. As a matter of fact, 
the main advantage of Eq.(15) is the clear separation between real and virtual 
contributions to the semihard cross section. More precisely, after summing, ac- 
cording with the AGK cutting rules, over all discontinuities of the semi-hard am- 
plitudes, which contribute to the inelastic process of interest, the product &i . . . &n 
is the remnant of the contribution from the real production terms. The product 
(1 — 6"tv + i) . . . (1 — ctq) is, on the contrary, the remnant of the contribution of the 
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virtual corrections[13]. The replacement &k — > OfeN m the former product, cor- 
responding to the real production process, allows one to generalize the functional 
in Eq.(ll) and to obtain the inclusive cross sections in the most general case of 
multiple parton interactions. One may therefore write 

e H [(3;z]= /EE^ W a ] «») "ViT } («i • • • O 

n m 

Q ( 16 ) 
x5 5] (S)®iN • • -®n{z](1 - &n+i) ... (1 - a Q )\[dudu' 

N=l ^ ' 



which gives the required inclusive cross sections via the relation 



z\ . . . dk n J 



d z (3- 



dk\...dk n J 5z(k\) . . .Sz(k n ) 
For later convenience 6h[/3;z] can also be expressed as 



(17) 



2 = 1 



6h N=feE } ^ • • • Un )~\ W< B l) «-p...u' m - 0) 

I Tl . TYl. 

n m 

n m 

x { n n i 1 + ®^ ^ ~ (is) 

n m 

-nn[ i -^ u *' u 5)]}n dudu ' 

i=i j=i 

We are now in a position to discuss the processes we are interested in, namely 
the events in which only two mini-jets are tagged. By setting n = 2 in Eq. (17) 
and using the second expression for @h[P', z], a lengthy but simple algebra yields 

9 



dk 1 dk 2 J dk x dk 2 

+ I D%\u,v)D%\u' -P,v'- P) d& ^ d& ^ dudu'dvdv> 
J 1 dk± dk 2 

+ f D^\u)D${u> -P,v>- d *(y') dudu > dv > 
J dki dk 2 

+ / D%\u,v)D%\u> - P) d ^^ d ^^ dudvdu> 
J dki dk 2 

where D^'(u) and D^ 2 \u, v ) are the one-body and two-body inclusive distribu- 
tions [10]: 



D^\u) =W^\u)+ j W^{u,u')du'+ i J W {i \u,u',u")du'du" + ... 

, -U2, u')du 

+ IJ w(4) ( u l> u ^ u> > u ") du ' du " ■■■ 

(20) 

In the r.h.s. of Eq. (19) every term has a clear physical interpretation. The 
first convolution is nothing but the usual single-collision contribution to the semi- 
hard cross section. The second term corresponds to two disconnected partonic 
collisions; finally, the last two entries correspond to those events in which a parton 
from hadron A or B has suffered a rescattering on hadron B or A respectively. 

From ref. [14] we know that the average number of rescatterings can be safely 
neglected in a typical hadron-hadron collision and for values of k m which allow 
the final state parton to be observed as an actual minijet in the final state. We are 
therefore allowed to neglect the last two terms in the r.h.s. of Eq. (19). The two- 
body inclusive distribution may be expressed by introducing the two body 
parton correlation C^: 
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D^( Ul ,u 2 ) = D^^D^M + ^C^( Ul ,u 2 ) (21) 

If one neglects both rescatterings and correlations in Eq.(19), one is left with the 
following simplified expression for the inclusive cross section: 



da% cl 
dk\dk 2 



d 2 f3 



where <g> is a compact notation for the convolutions appearing in Eq. (19). A 
possible further simplification follows from the assumption that D^'(u) has the 
factorized form 

D^(x,b) = f eff (x)F(b) (23) 
with the obvious normalizing condition 



J d 2 bF(b) 



1 (24) 
By substituting Eq. (23) in Eq. (22) one obtains 



daV} cl da s 1 da s da s 

H 37^37^ ( 25 ) 



dk\dk 2 dk\dk 2 a e ff dk\ dk 2 
where 

— = / d 2 (3 \ [ d 2 bF(b)F(b - (3)] 2 (26) 

(Jeff J U J 

and da s is the single collision expression, obtained by convoluting the elementary 
cross section with the usual one-body parton distribution f e ff(x). 

3. Numerical estimates 



11 



The formalism described in the previous section is a rather general approach 
to the problem of unitarity corrections in semi-hard interactions. Indeed the ex- 
pression for the inclusive cross section in Eq.(19) is completely general in the 
probabilistic picture of the semi-hard hadronic interaction. It is an exact conse- 
quence of the cross section as expressed in Eq.(14), which finds its justification in 
the AGK cutting rules [8]. In the inclusive cross section given by Eq.(19) all pos- 
sible multiple parton collisions are kept into account and multiparton correlations 
are treated at all orders. Consistently with the general principles, namely with 
the AGK cancellation [8, 15], the double inclusive cross section depends only on 
the single and double scattering terms. For a quantitative estimate of the role of 
unitarity corrections to the single scattering term, the required non perturbative 
input is represented both by the one-body parton distribution and by the 
two-body parton distribution D^ 2 \ The two-body parton distribution contains 
an independent information on the hadron structure with respect to actu- 
ally the two-body parton correlation C^ 2 \ While no experimental information is 
presently available on an indication is available from CDF on the scale fac- 
tor a e ff which characterizes the double parton interactions [16]. We will therefore 
limit our numerical analysis to the simplified case where is neglected and 
only disconnected parton collisions are taken into account, in such a way that the 
inclusive cross section is expressed by Eq.(25). All unitarity corrections to single 
scattering are therefore expressed by the second term in Eq.(25), which is obtained 
with the same input needed to evaluate the single scattering term, apart from the 
scale factor (J e ff, that summarizes all the geometrical details which enter in the 
unitarity correction. 

A few qualitative considerations are appropriate before illustrating the results 
of a quantitative analysis. By introducing the jet rapidities (y a , yt) and integrating 
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in the transverse momenta down to the lower cut off k m , while keeping fixed the 
azimuthal angle between the observed minijets 0, the inclusive cross section is 
expressed as: 

da 1 ^ 1 _ da s | 1 l_da^das_ 

d(pdy a dy b d(pdy a dy b a eff 2-rr dy a dy b 

In the limit of small relative rapidities y = y a — y b , a parton-parton interaction 
produces only two final state partons. Since they are back-to -back in 0, the single 
collision expression da s /d(j)dy a dy b is proportional to a Dirac delta 5((j) — %). This 
can be easily verified by setting y = in Eq.(7). On the opposite side, that is, 
for large values of y, the leading contribution to the r.h.s. of Eq.(7) comes from 
the n = term, for which the partons are decorrelated in (p. Physically, this 
is due to the large number of gluons radiated in the parton-parton interaction. 
Indeed, the flattening of the <f> distribution with increasing dijet rapidity gap was 
suggested [6] as a signature for the BFKL dynamics. From this point of view, a 
multiple partonic collision represents a background process which mimics the effect 
of multigluon emission. In the r.h.s. of Eq. (27), this background is described 
by the term weighted by the scale factor l/a e //. The experimental indication on 
the scale factor is 5.4 < a e ff < 29mb (90% C.L.)[16]. Unfortunately a e ff is not 
the only input variable which is still rather uncertain for a numerical computation. 
Indeed there is a large ambiguity already to compute the single scattering term. In 
fact to obtain the Lipatov's solution one needs to neglect the running of the strong 
coupling constant, in such a way that a s has to be considered as a parameter in 
the actual evaluation of a. Since the dependence of a on ct s , as it may be seen 
in Eq.(10), is exponential a numerical comparison of the two terms in Eq.(27) is 
rather uncertain. 
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To have a quantitiative feeling of the importance of the unitarity correction we 
have tried to estabilish a possible sensible choice of the input values of cr e // and 
a s by making a comparison with available experimental data. The experimental 
points in fig.l are the values of the cross section for production of minijets with 
k m > 5GeV measured by UA1[2]. The dashed curves refer to the single scattering 
integrated cross section with a s = .34 (upper curve) and a s = .29 (lower curve), 
corresponding to the values of the running coupling constant at the scale = 



k m / 3 and yCp = k m /2 respectively. The structure functions are the HMRS(B) 

structure functions [17]. The unitarized expression for the semihard cross section 
has a simple analytical representation when semi-hard rescatterings and multi- 
parton correlations are neglected[10]. Actually: 



where a s is the integrated single scattering inclusive cross section. The continuous 
curves in fig.l refer to the unitarized cross section an, as expressed in Eq.(28). For 
F(b) we have taken a gaussian, the width corresponding to a value of a e ff = 20mb. 
The two curves refer to the two different choices of a s mentioned above. The 
region identified by the two continuous lines contains the experimental points and 
therefore gives an indication on possible meaningful input parameters. One may 
also observe in fig.l how the rise of the experimental cross section is much closer 
to the rise of the unitarized curves than to the rise of the single scattering term 
alone. 

Before moving to different values of energy it is worthy to briefly comment on 
k m , which, to some extent, is a free parameter. A low value of k m corresponds to 
semi-hard cross sections that are well above a e ff (in the single collision approx- 
imation). In this conditions the contribution from multiple scatterings is largely 






(28) 
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dominant and the <p distribution is practically flat. On the contrary, large values 
of k m correspond to semi-hard cross sections that are negligible with respect to 
a e ff and no unitarity correction is required. Keeping this in mind, we realize that 
the interesting values of k m are those for which the total semi-hard cross section 
is comparable to cr e ff- This criterion yields km ~ 5.2 4- 6.1GeV at y/s = 1.8TeV 
and k m — 11.2 -4 12.7GeV at y/s = !8TeV, depending on the two different choices 
of values for a s which we have considered and for a e f / = 20mb. 

In order to have some quantitative indication on the effect that unitarization 
produces on the expectations based on the BFKL dynamics, we have studied the 
azimuthal correlation of the observed minijets, which, according with the BFKL 
dynamics has a distinctive dependence on the distance in rapidity. In fig. (2-a,b) 
we have plotted the differential cross section Eq.(27) as a function of </>, for fixed 
rapidities (y a ,yb) at y/s = 1.8TeV (a) and y/s = 18TeV (b) (the normalization 
is such that the curves take a value equal to unity at = and = 2ir). The 
naive 4> distribution, obtained by considering one elementary interaction only, is 
represented by the dashed line, while the continuous line describes the corrected 
distribution which takes into account an arbitrary number of parton-parton col- 
lisions. The flattening caused by the unitarity corrections is clearly visible: at 
y/s = 1.8TeV, fig. (2-a), the height of the central peak at = tt is reduced by a 
factor three approximately; the same trend, but with a stronger suppression of the 
correlation, is present at higher energies, see fig. (2-b). Figure (3) shows how the 
effect of unitarity corrections depends on the cutoff k m - By lowering this threshold 
we increase the semi- hard cross section and, accordingly, we enhance the probabil- 
ity of having several elementary parton collisions in each inelastic hadronic event. 
As a consequence, we expect the tagged minijets to become less and less corre- 
lated in the azimuthal angle <j>. This is confirmed by our plot which corresponds 
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to ^Js = 1.8TeV and to a rapidity gap y = 5, actually y a = 2.5 and y^ = —2.5. 
The different choices of the cut off k m are k m = 7GeV (solid line), k m = 6GeV 
(dashed line) and K m = bGeV (dotted line). It is worthwhile to stress that, for 
the lower choice k m = 5GeV, we cannot distinguish the distribution from a 
uniform one, unless we perform a quite accurate measure at the 3% level. Finally, 
fig. (4) shows how the correlation in the azimuthal angle of the tagging jets fades 
away as the rapidity interval is increased. 

4. Conclusions 

Minijet physics is the ideal tool to study BFKL dynamics. Indeed one comes 
closer and closer to the BFKL limiting case by keeping the lower threshold in 
transverse momentum k m of the observed minijets as small as possible. However 
the region of small k m is also the region where unitarity corrections become in- 
creasingly important. In the present paper we have made an attempt to estimate 
the unitarity corrections to the inclusive cross section for producing two minijets. 
After assuming the validity of the AGK cutting rules in semihard interactions, 
we have kept into account unitarity corrections by representing the hadronic pro- 
cess as a probabilistic superposition of multiple BFKL Pomeron exchanges. In 
the case of the inclusive cross section for producing two minijets, only the sin- 
gle and the double scattering terms contribute. With the purpose of making a 
quantitative estimate, we have considered the simplest possibility for the double 
scattering contribution. Actually we have neglected semi-hard parton rescatter- 
ings in the interaction, and two-body parton correlations in the two-body inclusive 
distributions. In this simplified case the unitarity correction depends on one sin- 
gle parameter only, namely cr e ff, that is the scale factor one needs to introduce 
in order to obtain the probability of the double interaction. For a quantitative 
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illustration of the effect of the correction term, a second parameter which has to 
be fixed is the strong coupling constant a s , whose value is not determined by the 
BFKL dynamics. Keeping into account the experimental suggestion on cr e // [16], 
we have fixed the input parameters by comparing with the UAl measurement of 
the semihard cross section for production of minijets[2]. Having selected in this 
way a possible range of values for the parameters, the indication we obtain from 
our numerical estimate is that at Tevatron energy the correction term to the single 
BFKL Pomeron exchange, depending on the actual quantity one is considering, 
may be larger than 100% for minijets with k m ~ 6GeV. When moving at LHC 
energies the same correction applies with values of k m ~ 12GeV. It is worthwhile 
noticing that, at Tevatron energy and with k m ~ 6GeV the average invariant mass 
of a partonic interaction is ~ .2TeV , while at LHC energies and with k m ~ 12GeV, 
the average invariant mass is ~ ITeV. The expectation is therefore that a sec- 
ondary BFKL Pomeron is exchanged in a large fraction of parton interactions at 
those values of invariant mass and at the corresponding hadron-hadron cm. en- 
ergy. A detailed experimental analysis of minijet production at Tevatron would 
therefore be of great importance both as a test of the BFKL approach, and to ac- 
cess the non perturbative information on the hadron structure which enters in the 
multiple parton interactions, whose detailed knowledge is of growing importance 
to understand hadron dynamics at higher energies. 
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Figure captions 

Fig. 1: Cross section for production of minijets with k m > 5GeV. Experimental data 
from UA1[2]. Dashed curves: single BFKL Pomeron exchange with a s — .34 
(upper curve) and a s = .29 (lower curve). Continuous curves: unitarized 
cross section, Eq.(28) in the text, same values of a s as in the previous case 
and a e ff = 20mb. 

Fig. 2: <j) distribution with unitarity corrections included (solid line) and in the single 
collision approximation (dashed line). N((f>) is proportional to the differential 
cross section given by eq. (27), with minijet rapidities kept fixed at y a = 2.5 
and yb = —2.5. The normalization is such that N(0) = N(2n) = 1. 

Fig. 3: (p distribution for several choices of the cutoff: k m i n = 7GeV (solid line), 
kmin = QGeV (dashed line), k m i n = 5GeV (dotted line). N(<f>) is defined as 
in Fig. (2) and unitarity corrections are included. 

Fig. 4: (j) distribution for different choices of the rapidity gap. The cutoff is k m i n = 
7GeV and the minijet rapidities are fixed at y a ^ = ±2.5 (solid line) and 
y a ^ = ±3.5 (dashed line). 
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